BARI-TH 444/02 
FIRENZE DFF 393/08/02 
UGVA-DPT-2002-09/1102 

Massive quark effects in two flavor color 
superconductors 

R. Casalbuoni°'^ F. De Fazio'^, R. Gatto^ G. Nardulli^'^ and 

M. Ruggieri^'^ 

Dipartimento di Fisica, Universita di Firenze, 1-50125 Firenze, Italia 
^I.N.F.N., Sezione di Firenze, 1-50125 Firenze, Italia 
Dipartimento di Fisica, Universita di Bari, 1-70124 Bari, Italia 
I.N.F.N., Sezione dz Ban, 1-70124 Ban, Italia 
Depart, de Physique Theorique, Universite de Geneve, CII-1211 Geneve 4, 

Suisse 

Abstract 

The high density effective theory formahsm (HDET) is employed 
to describe high density QCD with two massive flavors (2SC). The 
gap equation is derived and explicitly solved for the gap parame- 
ter. The parameters associated to the pseudo Nambu-Goldstone bo- 
son of U{1)a are evaluated in the limit /U ^ oo and m/fi fixed. 
In particular we find for the velocity of the NG boson the relation 

1 Introduction 



e 



In the present letter we study the effect of quark masses in the two flavor 
color superconductor (for a recent review see p[) in a BCS pairing. Cor- 
respondingly we will take the masses of two quarks not too different from 
each other. Massive quarks have been studied by several authors both for 
two and three flavors 0, 0, [Q, [Q. In most of the papers quark 
masses are treated perturbatively, under the hypothesis that they are much 
smaller than the energy gap. However this assumption is questionable for the 
strange quark, which may affect phenomenological analyses (for examples see 
Q, In the present paper we make a flrst attempt to go beyond the above 
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mentioned approximation by considering the limit — oo with x = m//i 
fixed. We will perform our analysis within the computational scheme of the 
High Density Effective Theory (HDET), |]TT] , |T2[ which is most suitable for 
calculations involving the high density limit of QCD. We will limit most of 
our analysis to the case of two flavors, though an extension to three flavors 
is in principle admissible. In Section 2 we establish our notations. Section 3 
is dedicated to the 2SC model with massive quarks where we study the gap 
equation. These results generalize the results of [0]. In Section 4 we analyze 
the properties of the pseudo Nambu Goldstone Boson (NGB) related to the 
spontaneous breaking of the U{1)a symmetry. We recover known results for 
the mass, although with a different method. On the other hand we find a 
new mass effect in the velocity of the pseudo NGB, which is given by 

2 v'/if -mfJ/ii - mi 

where fik,'^k [k = 1,2) are the masses and quark chemical potentials of the 
two quarks in the Cooper pair. This result is not in contradiction with the 
findings of other authors, see e.g. 0, since our corrections are of the order 
m/n rather than m/A. We also argue that this result should be true at all 
orders in our parameter expansion. 

2 The effective lagrangian for a massive gap- 
less quark 

Starting point is the QCD lagrangian 

^QCD = ^ (i -0 + /i 7° - ^ = ^ + /i 7° - (2) 
where we shall held 

x = — , 0<x<l (3) 
A'" 

fixed in the /i ^ oo limit. 

We introduce velocity dependent fields '^^{y) and '^^{y) 
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where A <^ 5 <^ /i and we have written the quark momentum as p'^ = 
fiv^ + i'^ , with v'^ = (0, v). In (^) we have used the projectors 

V± = . (5) 

These projectors can be obtained from the conventional ones, (1 ± a ■ v)/2, 
via a Cini-Touschek transformation [|T^. We recall that this transformations 
was introduced to study the ultra-relativistic limit of the Dirac equation. 
Owing to we can write 

Hy) = Y.(^"''"'[My) + '^v{y)] , (6) 

n 

where the sum represents an average over the velocity directions, with n = 
v/\v\ and arises because there are infinitely many ways to decompose the 
quark momentum according to (H). ipif{y) and "^viy) represent respectively 
the positive energy and negative energy solution of the Dirac equation asso- 
ciated to the lagrangian We now substitute in Cqcd and make use 
of the Riemann-Lebesgue lemma to eliminate one of the two sums appearing 
in the lagrangian (velocity superselection rule of HDET); therefore we get 

^o = T.{4^ 7°-^ + ^lit 7°^ + 2/i) + 1 7°-^ i^v + 4 ^ 7°^^.-} 

n 

(7) 

For quarks near the Fermi surface we can write |p| ~ fi\v\ and = //^ (1 — 
f^); therefore the parameter x defined in (|]) is also given by 

X = Vl -f2 . (8) 

We will use the identities 

7° = 7 -^ff = , (9) 

and 

V- a ■ vV+ = x{x - 7o)P+ . (10) 
We can integrate out in the functional integration \E'^ through the equation 

^^^= (11) 
iV ■ D ^ l[i 
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where we have used the notations = {l,v),V'^ = (1, —v) 

As in the case of massless flavor, the fields are of order with 
respect to ipij fields; therefore, at the leading order in one can write the 
Dirac part of the massive quark lagrangian as 

Co = J2 4^V-D^P.-. (12) 

n 

The integration of the antiquark fields gives a correction to the kinetic 
term of order This correction can be neglected in our limit. 

Since the sum over velocities is symmetric over the entire solid angle, one 
can write Cq as follows: 

^0= ^Y.{i^l{tV-D)i;+ + ijlitV-D)ij,} (13) 

n 

where the subscripts ± denote the velocity direction, i.e. 

V'± = i^±v ■ (14) 

This equation is formally the same as in the massless case, however one should 
remember that the velocity satisfies the condition (^. Another important 
difference is that we do not separate left and right degrees of freedom, and i/j± 
have two independent components. Finally the average over the directions is 
defined as follows: 

n 

The extra factor 1/2 in the average is introduced because, after the introduc- 
tion of the field with opposite velocity -0-, one doubles the degrees of freedom, 
which implies that the integration must be only over half solid angle. For 
the following it will be useful to introduce the Nambu-Gorkov notation 

3 2SC with massive quarks 

We consider the model with two massive quarks of different flavor in a color 
superconductor state 0. At the Fermi surface one has 

Pi = fJ'lVl, P2 = fJ'2V2 (17) 
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Let us put 

Hi = jj, — S/j, = fJ' + Sjj, . (18) 

The Cooper pair in the 2SC model couples quarks of different flavor, and the 
total momentum of the pair must vanish 

Pi+P2 = 0. (19) 

One gets immediately, from pi = p2 that 



^ Jl-xl + Jl- xl V ' 



where Xi = m^j jii. Eq. (|20| ) is a consequence of ([l9| ) and fixes 5// once the 
masses and yU are given. Here we are taking advantage of the fact that only 
strong interactions are considered and therefore has no constraints. By in- 
cluding also weak interactions, this line of reasoning might fail. As discussed 
in 1^, however, if the deviations from (|19|) and (pO|) are small, the value of 
the gap remains inaltered. Only for larger deviations one has a change, since 
there is a phase transition to the crystalline superconducting state |]14[, |jl5 



Eventually, for very large deviations, the normal non supercondcting state 
becomes energetically favorite. 

In order to write down the HDET lagrangian we first observe that the 
Dirac lagrangian not including the gap terms is simply given by two terms 
([l3|), one for each flavor. In the 2SC phase we have 

(0|V'^,C7 75V'ft.|0) = A6"'33e,,3 7^0, (21) 

where i and j are the flavor indices and a and (3 are color indices. Therefore 
we add the term 

^gap = -^A + -^A i (22) 

with 

>Ca = - E ^e^^' e.,3€,_ C 75^ft-,+ + h.c. (23) 

and 

n 

(24) 
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contains antiquark operators. These operators are formally suppressed 
but, as in the case of the Meissner mass |jl2[, they may give rise to divergent 
contributions. This happens for the mass of the U{1)a boson, where we 
get an unsuppressed term at the leading order in /x. We notice that by the 
equation we implicitly assume that the antiquark gap A is equal to the 



quark gap A. Using eqs. (11) and (|TOD we can write in momentum space for 
free quarks 

which is proportional to the quark mass. 

This introduces in the fermionic lagrangian in the 2SC model the non 
local interaction term 



£||Xi(xi + 7o) £||Xj(xj-7o) 

75 



- 75- 



l\\Xj{xj -7o) 



75- 



Vi\i2ni + Vi ■ 
i\\Xi{xi - 7o) 



Vi\{2fXi + Vi 



v,\{2^I, + Vri) 

+ h.c. (26) 



These terms will be taken in the large fi limit. Since can contribute, 
in the leading order in /i, only to divergent diagrams (i.e. giving rise to n 
factors), it may be treated as an insertion. 
Let us introduce 

x.= (j'-'' 

\Xa,2 . 

where Xa,i are the Nambu-Gorkov fields (see eq. 
Then we get: 



(27) 

TBI)) for the two flavors. 



(iVi-D'^'^ A"^ \ 

iVi ■ D*'"'^ A"^ 

-A"^ iV2-D'^^ 

V -A"^ Q iV2- Z^*'"^ ) 



Xfi 



A"^ = 75 A e"^^ . 
To get the quark propagator we define 

Di(£,A) = \/i-£\/2-£-A2, Z}2(^,A) = F2-£V1 



A^ 



(28) 
(29) 
(30) 
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We get 



(V2_ 








v.- 







D2 D2 







D2 





Do 











Vi 



(31) 



To get the gap equation in the 2SC model we use the Schwinger-Dyson 
equation with a four-fermion interaction: 



. 2 f 9fiuS' 



ab 



(27r)4 A2 



(32) 



where E = -[S''^ - Sp^], Sp^ stands for the inverse of the free quark prop- 
agator and 



\ 



i ■ Tf 



Tf being the usual 5'C/(3) generators. We get: 



—I 



A2 



V2-V^ (TaAT^T^ 47r(//'-5//') a(xi, X2) J 



\ 



iV2^ ■TJ'^PJ 



dH 1 
(27r)4i^i(£,A) 



(33) 



The factor (/i^ — (5/i^) a{xi,X2) = /'il^ip = /i2l'^2p = Pf implies a reduction 
of the phase space with respect to the massless case. For the sequel the 
following notations will be useful: 

a{xi,X2) ^ \Jl - xl^l - xl, l3{xi,X2) ^^l-xl-^l-xj. (34) 
The gap equation, for A 7^ is 



_ , r (r - 5^i^)a{xi, X2) V2 ■ Vi . 
' 3A%^ ^AA,a;i,a;2), 



(35) 
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where 



where the hmits of £\\ integration are ±5. One can easily check that for 
Xi = X2 = one recovers the gap equation for the 2SC model in the HDET 
approximation and for massless quarks[} Performing the integral we obtain: 



4a+^ S_ 

^'"^ [ilj^ -SiJ^)ag^ 2(1 + «) ) 

If the flavor 1 means the quark up and 2 the strange quark we can write 
Xi ^ 0; X = X2 and approximate (|20|) as follows: 



6fi 1 — vl — X 



or 



Ai 1 + Vl - a; 



(38) 



(39) 



and we get 



1-^ 7 r- (40) 

sinh 




In Fig. 1 we show the ratio of the gap for finite masses A over the gap 
for massless quarks Aq as a function of 6fi. In the plot X2 has been traded off 
for 6fi and different values of Xi have been used. Also the difference between 
the approximate expression and the exact one is given for xi = 0. 

The greatest effect of the strange quark mass on the condensation energy 
of the BCS pair is due to the reduction of the size of the Fermi surface. 
This implies that for small values of 6fi the condensation energy of the pair 
decreases linearly with 6fi. The xi = curve (dotted line) plots the same 
quantity as in ref. and agrees quantitatively with the results found there. 
It can be observed that also the results of 0, similarly to what we have 
done here, are neither obtained by expanding in m/A nor in m/fi. The 

iSee eqs. (2.188) and (2.184) of |l|. 
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10 20 30 40 

6^1 (MeV) 

Fig. 1 The ratio of the gap for finite masses A over the gap for massless 
quarks Aq as a function of The sohd hne refers to the approximate 
expression for A, see eq. (^OD, whereas the other three refer to the complete 
expression in eq. ( |37D with the dotted hne corresponding to Xi = 0, the 
dash-dotted hne to xi = 0.1, and the dashed one to xi = 0.25. The values 
used are /i = 400 MeV, A/g = 181 MeV, corresponding to Aq ^ 40 MeV. 
Notice that plotting 6fi is the same as plotting the squared masses difference, 
m| — ml due to the relation (|20|) . 

differences between our calculation and that of Rl follows. Our eq. 

(37) presents the advantage to offer an analytical expression; moreover it can 
be used when both masses, and not only one, are nonvanishing. On the other 
hand our results are only valid for rrii < fi; for higher values the evaluations 
of 0, that do not suffer of this limitation, should be used. 

We note that the insertion of would not change the main result we 
found, mainly the reduction of the phase space, because its effects would only 
change A by factors of order x/fi. 

4 U{1)a pseudo Nambu-Goldstone mode 

Besides color symmetry, that is broken from SU{3)c down to SU{2)c, the 
condensate (|2T|) also breaks U{1)b and U{l)em] however two linear combi- 
nations of these abelian subgroups remain unbroken and there are no NGBs 
associated to them. The vacuum expectation value ( pT]) breaks spontaneously 
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U{1)a, which is however also broken by the strong anomaly. At high /i this 
latter breaking is soft and one expects that the associated pseudo Goldstone 
boson, the 1]', is almost massless. In the sequel we will study the effect on rj' 
of the quark masses in the approximation of neglecting the strong anomaly. 

We introduce an external field a associated to the pseudo Goldstone boson 
and we put 



In ( p!6D we perform the substitution 

and Co+C/\ becomes Cq +>Ca+o-, with 

/ 

" J- n 

At the second order in the a field: 

U >^ l-\ Y 

and the quark-boson interaction lagrangian can be written as follows: 

J a fl Ja n 



(41) 



(42) 



(43) 



(44) 



where 



-1 



Vl 



1\ 



V-1 



1\ 



(45) 



The a kinetic term arises from loop expansion after functional integration of 
the fermionic degrees of freedom. Therefore we have the following effective 
action 



1 



"^eff = i-Tr dxdy <iS{y,x 



2ia{x] 



AiToiS{x,y) 



fa 
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+ zTr J dx l^iS{x,x)^^Aa^{x)tTij . (46) 

In evaluating the trace over the spin indices one has to include the V+ pro- 
jector which gives a factor 2 if no other 7-matrices are involved. The two 
terms correspond to the self-energy and a tadpole diagram respectively, see 
fig. 2. 



Fig. 2 One-loop diagrams. External lines represent the pseudo Goldstone 
field. Full lines are fermion propagators. 

The result of the calculation of the effective lagrangian in momentum 
space is as follows: 



16/ii/i2aA^ ^ f ,2 

J a n 



h 1 ^ 



Di(£) Di(£ + p) 



(47) 



tad 



7r3/2 ^ } \D,{e,A) D2(^,A), 
We find 

Cs.e.iP = 0) + ^tad = . (49) 

We have to add the contribution from the antiquark gaps. The pseudo 
NGB masses in color superconducting phases take contribution from different 
sources. In the case of the singlet pseudo NGB, as shown in 0, its mass arises 
from an effective four quark operator obtained by integrating out the electric 
gluon field; the result is therefore proportional^ to mim2A^//i^. In this 
approach the antiquark gap terms do not appear explicitly as the antiquark 
field operators are integrated out. 

^The calculations in ||] are in the Color Flavor Locking phase and should be adapted 
to the 2SC case. 



11 



Also in the present approach we integrate out the antiquark operators, 
as discussed above and is expressed in terms of the quark fields only. It 
contributes both to £s.e.(0) and to Ctad, through two additional vertices axx 
and aaxx that can be constructed analogously to using (p6D and (^l]). 
We insert these vertices in the expressions for the self-energy and for the 
tadpole once, since we will take only the leading term in the quark masses 
to the mass of the Goldstone. The two contributions do not cancel out and 
one finds 



^>^e//(0) = — ^7^X1X2/ . (50) 



For / we get: 

dH J^im « ^2..,^ log + . 



This is the leading contribution and arises from the first term in (p6D, since 
the other two terms are of order /i^ and not /i^ log /i. 
Let us then consider £e//(p) — 'Ce//(0). 



{[{V, ■ p)af [V, ■ if + \iy, ■ v)of [y^ (1 ^ 2) . (51) 



One can show that the second term on the r.h.s. of ( pi]) vanishes, while one 
finds 

(52) 



We get therefore 



4 /ii/i2 a 



TT (4a + p^) ' ,j 
From which one finally gets, in coordinate space 
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where to get the canonical normahzation for the kinetic lagrangian one has 
to put 

2 _ 8/ii/i2 \Vl\\v2\ f^^. 

Using the resuh in (|50|) one finds at the leading order: 

m2=4A2 ^log-^, (54) 

which coincides with the result of 0. It should be noted that in the 
masses of the NGB bosons in the CFL case do not present any logarithmic 
enhancement factor. The major difference between our approach and that 
of IP], a part from the difference in the models that are studied (CFL ver- 
sus 2SC), is in the coupling between quarks, as we use a BCS four-fermion 
interaction while in |§] one considers one gluon exchange interaction. 

Nonzero quark masses produce a change in the coupling constant and 
a reduction of the velocity of the pseudo-NGB in the medium, i.e. 



1 — — x 



^2 



whereas for massless quarks v"^ = 1/3. We wish to stress that Eq. 
constitutes a new result which goes beyond the expansion in m/A commonly 
used in the computation of mass effects in color superconductivity. 
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